The PAMELA and ATIC Excesses From a Nearby Clump of Neutralino Dark Matter 
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In this letter, we suggest that a nearby clump of 600-1000 GeV neutralinos may be responsible 
for the excesses recently observed in the cosmic ray positron and electron spectra by the PAMELA 
and ATIC experiments. Although neutralino dark matter annihilating throughout the halo of the 
Milky Way is predicted to produce a softer spectrum than is observed, and violate constraints from 
cosmic ray antiproton measurements, a large nearby (within 1-2 kiloparsecs of the Solar System) 
clump of annihilating neutralinos can lead to a spectrum which is consistent with PAMELA and 
ATIC, while also producing an acceptable antiproton flux. Furthermore, the presence of a large dark 
matter clump can potentially accommodate the very large annihilation rate required to produce the 
PAMELA and ATIC signals. 
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Recent observations by the PAMELA pLj] and ATIC ^ 
experiments have revealed a surprisingly large flux of 
high energy electrons and positrons in the cosmic ray 
spectrum. These observations strongly indicate the pres- 
ence of a relatively local source of energetic pairs. Al- 
though the origin of these particles remains unknown, a 
nearby pulsar [S, A] and dark matter annihilations 0, @, 
0, M, S [ill. Il2j have each been proposed as possible 
sources. 

Efforts to explain the PAMELA and ATIC excesses 
vifith annihilating dark matter have faced a number of 
challenges, however. In particular: 

• The spectrum of electrons and positrons predicted 
to be generated in the annihilations of most dark 
matter candidates is much too soft to fit the ob- 
servations of PAMELA and ATIC. If WIMPs an- 
nihilating throughout the halo of the Milky Way 
are to produce the spectral shape observed by 
these experiments, they must annihilate mostly to 
charged leptons {e~^e~, jX^ fi~ ^ and/or r+r") Q. 
While models have been proposed in which this 
is the case 0, many of the most often studied 
WIMP candidates (including neutralinos) are pre- 
dicted to annihilate dominantly to quarks and/or 
gauge bosons [isj. 

• The dark matter annihilation rate that is required 
to generate the observed spectrum of cosmic ray 
electrons and positrons is much higher (by a factor 
of ^lO^-lO"^) than is predicted for a typical thermal 
relic distributed smoothly throughout the Galactic 
Halo [3]. To normalize the annihilation rate to the 
PAMELA and ATIC signals, we must require ei- 
ther large inhomogeneities in the dark matter dis- 
tribution which lead to a considerably enhanced 
annihilation rate (ie. a " boost factor"), and/or 
dark matter particles which possess a considerably 
larger annihilation cross section than is required 



of a thermal relic. This latter option requires ei- 
ther a non-thermal production mechanism in the 
early universe, or an enhancement of the annihila- 
tion cross section at low velocities, such as through 
the Sommerfeld effect 0,[i|. 

• The very large annihilation rate required to gen- 
erate the PAMELA and ATIC signals can also 
lead to the overproduction of cosmic ray antipro- 
l ld . llj . gamma rays, and synchrotron emis- 




For example, Ref. [ll] finds that a 1 TeV 
WIMP annihilating to W^W~ would be expected 
to exceed the observed cosmic ray antiproton flux 
by a factor of approximately 5-10 if the overall 
annihilation rate is normalized to the PAMELA 
positron fraction. 

Collectively, these considerations appear to strongly 
limit the type of dark matter particle that could poten- 
tially generate the observed spectra of cosmic ray elec- 
trons and positrons. Under the astrophysical assump- 
tions that are typically adopted, dark matter must anni- 
hilate almost entirely to charged leptons, and at a very 
high rate, if they are to accommodate the observations 
of PAMELA and ATIC. 

Here, we propose an alternative scenario that can alle- 
viate all three of the challenges listed above. In particu- 
lar, we consider the case in which the high energy cosmic 
ray electron and positron spectra are dominated by neu- 
tralino annihilations taking place in a nearby clump of 
dark matter. In such a scenario, WIMPs which anni- 
hilate to non-leptonic final states can still generate the 
spectra observed by PAMELA and ATIC, without vio- 
lating constraints from cosmic ray antiprotons, gamma 
rays, or synchrotron emission. 

As electrons and positrons propagate through the radi- 
ation flelds and magnetic field of the Milky Way, they lose 
energy through inverse Compton and synchrotron pro- 
cesses. As a result, the spectrum from distant objects 
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is softened relative to that originating from more local 
sources. This is illustrated in the top frame of Fig. [U 
where we show the spectrum of electrons plus positrons 
from the annihilations of WIMPs in a stationary clump 
0.1, 1, 2, or 4 kpc from the Solar System. Here, we 
have considered an 800 GeV WIMP which annihilates 
to W'^W~ (such as a wino-like neutralino, for example). 
We compare this to the spectral shape predicted from an- 
nihilations throughout a smooth (Navarro-Frenk- White, 
NFW [l^) halo profile, and to the spectral shape of elec- 
trons/positrons prior to propagation. From this, it is 
clear that a nearby clump of annihilating dark matter 
can lead to a spectrum that is considerably harder than 
is predicted from annihilations throughout the halo at 
large. 

To calculate the cosmic ray electron/positron spectrum 
taking into account the effects of diffusion and energy 
losses, we have solved the propagation equation [15|: 
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where dNe/dEg is the number density of positrons per 
unit energy, K{Ee,x) is the diffusion coefficient, and 



h{Ee 



is the rate of energy loss. The source term. 



Q{Ee,x), reflects both the distribution of dark matter 
in the Galaxy, and the mass, annihilation cross section, 
and dominant annihilation channels of the neutralino. 
Throughout this letter, we adopt the following diffusion 
parameters: K{Ec) = 5.3 x lO^^ (£;c/4 GeV)°"3 cm^/s, 
and b{Ee) = IQ-i^ (E^e/l GeV)^ s'^ We also se- 
lect boundary conditions corresponding to a slab of 
half-thickness 4 kpc, beyond which cosmic ray elec- 
trons/positrons are allowed to freely escape the Galactic 
Magnetic Field. 

Dark matter substructures close to the center of a dark 
matter halo, such as near the Solar System in the Milky 
Way, typically move with velocities of a few times the ro- 
tational velocity (wrot = 220km/s for the Sun). We can 
expect a broad distribution of clump velocities, but typ- 
ical values are approximately 400km/s. The stationary 
solution to Eq. ([f]), such as used in Refs. [H, [13], is a 
good approximation when vR ^ K[Ec\, where R is the 
distance to the clump and v is the velocity of the clump 
with respect to the ISM. Given our choice, K{Ee) « 
685 (i?o/100 GeV)"-^^ kpckm/s, we see that clump mo- 
tions are largely unimportant for (i?c/100 GeV)°'^^ S> 
i?/kpc, but increasingly important at lower energies. In 
Fig. 1 we show the electron plus positron distribution 
from dark matter annihilations for a stationary clump 
and a moving one. Including the motion of the clump 
hardens the spectrum further. At energies below a few 
hundred GeV, the spectrum is suppressed relative to the 
steady state case. In contrast, the spectrum is unchanged 
above a few hundred GeV. 
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FIG. 1: The electron plus positron spectrum from a clump 
of 800 GeV neutralinos annihilating to W^W' . In the top 
frame, the clump is assumed to be stationary and at distances 
of 0.1, 1, 2 or 4 kpc from the Solar System. For comparison, 
we show with arbitrary normalization the spectral shape prior 
to propagation (dotted line) and the result for a smooth NFW 
halo profile (dashed line). In the lower frame, results are 
shown for a clump moving at 400 km/s relative to the Solar 
System. The thick solid line denote the case in which the 
clump has recently reached the Solar System. Other line types 
describe the case in which a clump which has passed through 
the Solar System and is now 1 or 2 kpc away (dotted, top-to- 
bottom), a clump which is approaching the Solar System and 
is currently 1 or 2 kpc away (dashed, top-to-bottom), and a 
clump which passed with a closest approach of 1 kpc away 
from the Solar System and is now V2 or 2^2 kpc away (thin 
solid). In each case shown, an annihilation rate of 2 x 10'^^ 
s~^ was used. 



In Fig. [21 we compare our results to the measure- 
ments of PAMELA (top) and ATIC (bottom). Again, 
we have considered a dark matter particle annihilating 
to W^W~, such as a wino-like neutralino. From these 
frames, it is clear that a nearby clump of such particles 
can accommodate the measurements reported by these 
experiments. As the ATIC feature becomes less promi- 
nent for more distant clumps, the origin of the signal 
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FIG. 2: The positron fraction (top) and the electron plus 
positron spectrum (bottom) from a nearby clump of anni- 
hilating neutralino dark matter. The solid lines denote the 
result for a very nearby clump of 600 GeV neutralinos, while 
the dashed and dotted lines correspond to 800 GeV neutrali- 
nos in a clump 1 and 1.2 kpc away, respectively. To normalize 
the curves, we used approximate annihilation rates of 7 x 10"^^ , 
1.7 X 10^'', and 2.2 x 10^^ per second for the D = 0, 1, and 1.2 
cases, respectively. Each case provides a good fit to both the 
PAMELA and ATIC data. For comparison, we show in the 
top frame as a dot-dashed line the astrophysical expectation 
from cosmic ray secondary production \18^ . 



must lie within approximately 1 or 2 kpc of the Solar 
System to generate the observed spectral shape. 

To normalize the signals shown in Fig. ^ we have used 
annihilation rates of 7 xl03^ 1.7x10^^, and2.2xl03'^ per 
second for clumps at distances of 0, 1 and 1.2 kpc, respec- 
tively. These very large annihilation rates require us to 
consider a very large or dense clump of dark matter. For 
comparison, consider the dwarf spheroidal galaxy Draco. 
Assuming a smooth NFW halo, dark matter particles 
could be annihilating within Draco at a rate as large as 



5 x 1035 s-i (600 GeV/TOx)^ (o-«/3 x 10" 



mately 3 x 10^ solar masses, a Draco-like object located 
within a fraction of a kpc from the Solar System would be 
only marginally capable of generating the PAMELA and 
ATIC signals for a thermal WIMP distributed through- 
out the clump with a smooth NFW profile. 

The mass and density required of the clump could be 
relaxed considerably, however, if the neutralino were pro- 
duced non-thermally in the early universe. In particular, 
an 600 GeV wino-like neutralino is predicted to have an 
annihilation cross section of cru « 3 x 10"^^ cm'^/s, which 
is approximately an order of magnitude larger than the 
value predicted for a typical thermal relic. Alternatively, 
inhomogeneities within the clump's dark matter distribu- 
tion can lead to a higher annihilation rate. If some com- 
bination of these effects enhance the annihilation rate by 
a factor of '--^lO, then the calculations of Ref. [l^ can be 
used to estimate that there is a 0.2% to 0.02% chance 
of a sufhciently massive and/or dense clump being close 
enough to the Solar System to generate the observed cos- 
mic ray electrons and positrons. 

Lastly, we turn our attention to the constraints that 
can be placed from observations of gamma rays, syn- 
chrotron emission, and cosmic ray antiprotons. For a 
typical distribution of dark matter throughout the halo 
of the Milky Way, the strongest indirect detection con- 
straints come from limits on the gamma ray [20| and 
synchrotron 2l| flux from the inner galaxy. In the case 
in which the PAMELA and ATIC signals are dominated 
from nearby annihilations, the annihilation rate in the 
inner Milky Way can easily be low enough to evade these 
constraints. A nearby clump of neutralinos would, how- 
ever, produce an extended source of gamma rays which 
could potentially have been observed by EGRET (or in 
the future by FERMI/GLAST). Considering a clump of 
800 GeV neutralinos annihilating at a rate of 1.7 x 10"^^ 
s~^, 1 kpc from the Solar System, we estimate the flux 



of gamma rays at 1 GeV to be E^dN^/dE. 
10-6 GeVcm-2 



5 X 



-sr i(200pc/r)2, 
x 10-5 GeV cm-2 



and at 10 GeV to be 



Thus, even with its very substantial mass of approxi- 



E?^dN^/dE^ « 2 x 10-5 GeV cm-2 s-i sr-i (200pc/r)2, 
where r is the radial size of the annihilating region. As 
this is below or comparable to the diffuse background 
measured by EGRET , the presence of such a clump 
appears to be consistent with this constraint. 

In addition to measuring the positron fraction, 
PAMELA has also published their measurement of 
the antiproton-to-proton ratio in the cosmic ray spec- 
trum [23|. This result, which is consistent with pure 
secondary production of antiprotons, can be used to con- 
strain the rate of dark matter annihilations in the Milky 
Way halo. Ref. [Il|, for example, found that if TeV-mass 
WIMPs annihilating to W^W~ throughout the halo of 
the Milky Way is normalized to produce the PAMELA 
positron fraction, then the flux of antiprotons would ex- 
ceed the observed antiproton-to-proton ratio by about an 
order of magnitude. This conclusion is altered consider- 
ably in the case of a nearby annihilating clump, how- 
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ever. 600 GeV (100 GeV) electrons and positrons typ- 
ically lose most of their energy by the time they have 
traveled ~ 1 kpc (a few kpc). In contrast, antiprotons 
propagate without significant energy losses, and can con- 
tribute to the cosmic ray spectrum even if they originate 
from very distant regions of the galaxy. As a result, the 
constraint from cosmic ray antiprotons is relaxed consid- 
erably when only the local annihilation rate is normal- 
ized to the PAMELA positron fraction. Furthermore, 
one should keep in mind that modifications to the propa- 
gation model (especially the width of the diffusion zone) 
could lead to an antiproton constraint which is less strin- 
gent by a factor of 20 or more than described in Ref. [lH . 

In summary, we have discussed in this letter the pos- 
sibility that the excesses recently observed by PAMELA 
and ATIC may be the product of a nearby clump of an- 
nihilating neutralino dark matter. This scenario solves 
three problems that are typically faced when attempt- 
ing to explain these signals with annihilating neutralinos. 
First, the spectrum is hardened considerably (especially 
when the motion of the clump is properly accounted for), 
bringing the predicted spectrum into line with the mea- 
surements of PAMELA and ATIC. Second, a larger an- 
nihilation rate is expected, making the observed flux less 
challenging to accomodate. Third, the constraints from 
cosmic ray antiprotons are relaxed considerably, as well 
as those from gamma rays and synchrotron emission from 
the inner Milky Way. 

To determine whether or not a nearby clump of anni- 
hilating neutralinos is in fact responsible for the ATIC 
and PAMELA signals, further data will be required. In 
addition to more data from PAMELA (at higher ener- 
gies, and with greater exposure), ground based gamma 
ray telescopes such as HESS and VERITAS should be 
capable of measuring the cosmic ray electron spectrum 
over the energy range of the ATIC feature with higher 
precision than is currently available (23 | . 

Data from the FERMl/GLAST gamma ray telescope 
will be valuable in testing this hypothesis in two ways. 
Firstly, this experiment may be capable of detecting 
gamma rays from the clump itself, which would provide a 
smoking gun for the scenario described in this letter. Sec- 
ondly, FERMI will also detect very large nubmers of cos- 
mic ray electrons/positrons, which could potentially en- 
able the detection of a small (~0.1%) dipole anisotropy in 
their angular distribution 0]. Such an anisotropy would 
not be expected if dark matter annihilating throughout 
the halo were responsible for the PAMELA and ATIC 
signals. 
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